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In this paper, we derive various integral equations related to the 
scattering of time-harmonic electromagnetic fields by perfect conduc- 
tors with 2-dimensional geometry. The fields may be expressed in 
terms of two solutions of a scalar wave equation and decomposed 
into E waves and H waves. We consider the case in which part, or 
all, of each of the perfectly conducting cylindrical scatterers may be 
infinitesimally thin, and show that a standard integral equation, 
used in the case of H waves, does not determine the current density 
on the infinitesimally thin parts of the scatterers. We derive an 
alternate integral equation which does not suffer from this defect. 
This equation has been used by J. L. Blue in the numerical solution 
of the problem of scattering by an infinitesimally thin strip. 

I. INTRODUCTION 

In this paper, we consider the scattering of time-harmonic electro- 
magnetic fields by perfect conductors with 2-dimensional geometry, in 
which the boundaries are independent of the 2-coordinate. The z- 
dependence of the fields is assumed to be of the form exp(ik sin a z), 
where k is the free space wave number and | a \ < tt/2, so that the 
scattering of obliquely incident plane waves may be investigated. It is 
known 1 that the electromagnetic fields may be expressed in terms of 
the longitudinal components, E z and H z , and that each of these two 
quantities satisfies the scalar wave equation, with wave number 
k cos a. Moreover, since the boundary conditions on a perfectly 
conducting surface imply that both E z and the normal derivative of H z 
are zero, there is no coupling between E z and H z , and we refer to E 
waves and H waves, respectively. 

Integral equations for scattering problems have been considered by 
numerous authors. A relatively recent treatment of this topic is that of 
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Poggio and Miller, 2 but they give only a brief discussion of the 2- 
dimensional case. A useful discussion of integral equations for the 
scalar problem is given by Noble. 3 Poggio and Miller state that the 
integral equation which they derive for the surface current in the case 
of H waves is useless when the scatterer is infinitely thin. Noble points 
out that the corresponding integral equation, when applied to the 
problem of scattering by an elliptic cylinder, degenerates as the eccen- 
tricity tends to unity, so that the scatterer becomes an infinitesimally 
thin strip. 

In this paper, we consider the case in which part, or all, of each of 
the perfectly conducting cylindrical scatterers may be infinitely thin. 
We derive an alternate integral equation for the current density on the 
scatterers, in the case of H waves, which does not degenerate on the 
infinitesimally thin segments. We discuss the relationship between this 
integral equation and one derived by Mitzner. 4 Our integral equation 
has been used by Blue 6 in the numerical solution of the problem of 
scattering by an infinitesimally thin strip. We also point out how the 
integral equation which does degenerate may be used to calculate H z 
on both sides of the infinitesimally thin segments, once the entire 
current density on the scatterers is known. An integral equation which 
degenerates in the case of E waves is also derived, and this may be 
used analogously to calculate the values of the normal derivatives of 
E z on both sides of the infinitesimally thin segments, once the entire 
current density on the scatterers is known. 

In Section II we briefly derive expressions for the transverse com- 
ponents of the field in terms of the longitudinal components E, and 
H z , and show that the latter quantities both satisfy a scalar wave 
equation. We also derive the boundary conditions on a perfectly 
conducting surface, and give an expression for the current density on 
the surface. The total fields are expressed as the sum of the incident 
and scattered fields. In Section HI we derive an integral representation 
for the scattered field in terms of the total field, and its normal 
derivative, on the scatterers. A nondegenerate integral equation for 
the current density on the scatterers is obtained in the case of E waves, 
by using this representation as a point on the boundary is approached. 
In the case of H waves, it is shown that the corresponding integral 
equation degenerates on the infinitesimally thin segments, since it 
contains an unknown quantity besides the current density. 

In Section IV we derive representations for the transverse compo- 
nent of the gradient of the scattered field. By using this representation 
to calculate the normal derivative of the scattered field as a point on 
the boundary is approached, we obtain a nondegenerate integral 
equation for the current density on the scatterers in the case of H 
waves. In the case of E waves, it is shown that the corresponding 
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integral equation degenerates on the infinitesimally thin segments. 
The implications of these results are discussed. 

II. THE ELECTROMAGNETIC FIELDS 

We first write down equations which describe the electromagnetic 
fields due to scattering by perfect conductors with 2-dimensional 
geometry. If we suppress the factor exp(— iut), where w is the angular 
frequency, the divergenceless electric and magnetic fields E and H, in 
free space, satisfy Maxwell's equations 6 

V X E = iwjUoH, V X H = -iweoE, (1) 

where juo is the permeability and eo is the dielectric constant. The free 
space wave number is k = w(/ioeo) 1/2 . We consider the case of a 2- 
dimensional geometry in which the boundaries are independent of the 
coordinate z, and assume that the 2-dependence of the fields is of the 
form exp(ik sin a z), where | a | < it/2. This will allow us to consider 
scattering of obliquely incident plane waves. Accordingly, we now 
suppress the factor exp(ik sin a z), and write 

V = V t + ikainaiz, E = E, + E z \ z , H = H, + H z i z , (2) 

where i z is a unit vector in the z-direction, and the subscript t refers to 
the transverse components. 

If we split eqs. (1) into their transverse and longitudinal components, 
we obtain 

V,E Z X i z - ik sin a E, X i 2 + wo/toH,, (3) 

V«ft X \ z - ik sin a H, X i 2 - iwcoE,, (4) 

and 

V,xE,= iu3HoH z i z , V,xH, = -iw€ E z i z . (5) 

It is convenient to define the transverse wave number k t = k cos a. 
Then, from (3) and (4), it follows that 

ktEt = i(k sin a V t E* + wpovYH, X i 2 ) (6) 

and 

kfHt = i(k sin a V,H Z - ue V t E z X i z ). (7) 

Hence the transverse fields are expressed in terms of the longitudinal 
components. If we substitute these expressions for E, and H ( into (5), 
we obtain 

(V? + #)E 2 = 0, (V? + A?)H z = 0, (8) 

where we have used the relationships 7 

V,x(V f V)-0, V, x(V,Vxi 2 ) = -V?Vi z . (9) 

INTEGRAL EQUATIONS FOR SCATTERING PROBLEMS 411 



Hence, as is known, 1 the longitudinal components of the field satisfy 
the scalar reduced wave equation. 

The boundary conditions 8 on a perfectly conducting surface are that 
the tangential components of the electric field and the normal com- 
ponent of the magnetic field vanish, i.e., 

Exn = 0, Hn = 0, (10) 

where n is a unit vector normal to the surface, directed into the 
scattering region. Because of the 2-dimensional geometry, n-i 2 = 0. 
From (2), (6) and (7), it follows that the boundary conditions (10) are 
equivalent to 

flH 
E z = 0, — - = on a perfectly conducting surface. (11) 
dn 

The current density on the surface 8 is K = n X H. We assume that t, 
n and i 2 form a right-handed set of unit vectors. Then, from (2) and 
(7) we find that 

if BE Z . . dH : 



kt \ dn ds 

where s denotes arc length along the cross-sectional boundary curve, 
and t is a unit vector tangent to the curve. 
We write the total fields as the sum of incident and scattered fields, 

E = E' + E 8 , H = H' + H 8 . (13) 

In the case of incident plane waves we have, in Cartesian coordinates 
(x, y, z), 

E z = E exp[ik t (x cos + y sin /?)], H z = 0, (14) 

for an E wave, and 

H z = Ho exp[ik t (x co&fi + y sin p)], E l z = 0, (15) 

for an H wave. The factor exp[i(k sin a z - ut)] has been suppressed. 

III. INTEGRAL EQUATIONS DERIVED FROM REPRESENTATION FOR 
THE SCATTERED FIELD 

We first derive some integral representations for the scattered fields. 
We suppose that 

(V? + A:?W = 0, (V? + khr = 0, (16) 

and set i//' = E z and \p 8 = E\, or ^' = H z and \p" = HI, corresponding 
to E waves, or H waves, respectively. We let 

*-** + *•■ (17) 
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and then the corresponding boundary conditions are, from (11), 

dxP 

$ = 0, or — = 0, on a perfectly conducting boundary. (18) 
an 

We introduce the 2-dimensional Green's function 3 

G{T,p)=- A m ) (k t R), (19) 

4 

where H§ ] denotes a Hankel function 9 of zero order, and 

r = xi x + yi y , p = £i x + v i y , R = r - p, i? = |R|, (20) 

where i x and i^ are unit vectors in the x- and .y-directions. Then 3 

(V? + #)G = 0, R*0, (21) 

We consider the case in which part, or all, of each perfectly con- 
ducting cylinder may be infinitesimally thin, and let the cross-sectional 
boundary curve of the jth cylinder be denoted by Cj = Tj U Lj U Lj, 
where Lf and Lj denote opposite sides of the infinitesimally thin 
segment(s) Lj. The segments Lj may be disjoint, as may be Tj also, as 
depicted in Fig. 1. The curves Cj are assumed to be piece wise differ- 
entiable. We consider a point r exterior to all the curves Cj, and apply 
Green's theorem 10 in the region A exterior to the curves Cj, exterior to 
| p — r | = e, and interior to | p | = t, as depicted in Fig. 1. Then, with 
C = UjCj, 



>Cu\ P -r\-tU\ P \-T - 



f [ BG BV ] 

*>) — -G^-(p) 

Jculo-H-.U.nl-, - dH dfl J 

1 



r(p) — -G-^( P )\ds 



= (yViG-GVfy 8 )dA = 0, (22) 



from (16) and (21). Because of our choice of n, the normal derivatives 
are directed into the region A. 
Now, from (19), since 9 

Htf } (k t R) = — \og(k t R) + 0(1), for k t R « 1, (23) 

IT 



it follows that 



lim I [^ 8 ^- G ^' d; =^ ir] {24} 

\p-r\-f 



Also, since 11 

/ 2 \ 1/2 / r w "i\ 

for k,R»l, (25) 



HltHk.R)- (~y exp(z 



k t R-\ 
4 
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Fig. 1 — Cross section of cylindrical scatterers. 



BnJ LsdG_ G df\ ds= (26) 



if the scattered fields satisfy the radiation condition 12 



limp 1/2 (^-^ 8 H0> ( P = |p|), 

p_»oo \ dp 



(27) 



which we assume to be the case. If we let € — > and t — ► oo in (22), it 
follows from (24) and (26) that 

^(r) = J[^(p)^-G^(p)]rf S . (28) 

If we consider the incident field \p\ and apply Green's theorem in 
the region(s) enclosed by Tj, and use (16) and (21), we obtain 
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Also, because of the continuity of \p' and G, and their normal deriva- 
tives, on Lj, it follows that 



J J 



dG d^ 

on on 



ds = 0, 



(30) 



since the normals are reversed on opposite sides of Lj. Hence, since 

c=u,(r,uz,;uL7), 



*'<»£-<4«» 



ds = 0. 



(31) 



Consequently, with the help of (17), we may rewrite (28) in the form 



^ 8 (r) = 



c - 



dG „ d\b 



ds. 



(32) 



The advantage of doing this is that, because of the boundary conditions 
in (18), one of the two terms in the integrand in (32) vanishes. 

We first consider the case of E waves. Then, from (11), (19), and 
(32), we have 



dE 2 



El(r) = -7 HP(ktR) - 1 (p) ds. 



dn 



If we now let r tend to a point on C, we obtain 



E\ 



dE z 



) c dn 



«»-JJ[ 

since E* z + E\ = on C. We may rewrite (34) in the form 

i f dF 

Ei(r)=- HS*ihR)-j£(p)da 

BE 



(33) 



(34) 



+ l - m i) (k t R) 



dn 



(P) 



BE, 

■7— (P) 

dn 



ds, reC, (35) 



where T = Uyly and L = U/L>. This integral equation may be used to 
determine the current density on T and the total current density on L, 
which suffices to determine the scattered field from (33). However, 
(35) does not yield the separate values of the normal derivative of E z 
on either side of L. We will return to this point in the next section. 
We now consider the case of H waves. Since 9 



— Htf ) (k l R) = -k l H[ 1) (k t R), 
dK 



(36) 
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it follows from (11), (19), (20), and (32) that 



Hl(r) = l -k t j 



H[ 1) (k l R)?^H z ip)ds. (37) 



Let ro be a point on T, not at a corner, and let a be a small segment of 
T containing r . Then, as seen from Fig. 2, by considering the angle 
&f> subtended by the element 8s of a, and letting 8s — * 0, we obtain 

Since, from (23) and (36), k,RH[ l) (k t R) -» -2i/m as k,R -> 0, it follows 
that 

ton ±k t j HFikB) ^ H Z (P) ds ~ ^- j d*. (39) 

But since we have assumed that ro is not at a corner, the angle 
subtended at ro by a tends to tt as the length of a tends to zero. Hence, 
from (37), we have 



(40) 



Hl(r) = i H 2 (r) + £ k t P J H[ n (ktR) ^ ft(p) ds, r e V, 

where P denotes the principal value of the integral, corresponding to 
the limit of the integral over C — a as the length of a tends to zero, and 
T' denotes T less its corners. Since H z = H l z + H%, we may rewrite (40) 
in the form 

\ H z (r) = HUr) + - d k t P \ H[ 1] (k t R) — H z (p) ds 
2 4 J 7i 



i*I 



£^(ik»A) I . {n + [^ 2 (p)] + 



+ n_[if z (p)]_} rfs, reT'. (41) 

Now let r be a point on Lj, not at a corner (or endpoint), and let a 
be a small segment of Lj containing ro. Then, from (38), if we let r 
tend to r from the Lj side, we obtain 

lim -k t j Hi'HktR)^- {n + [H z (p)] + +n-[H z {p)]-} ds 

.[r ] + 4 J a K 

-jttftWJf-WWWJ*. (42) 



r-»i 



since n_ = — n+. Hence from (37), it follows that 
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Fig. 2 — Angle subtended by an element of a cross-sectional boundary curve. 

l -k t J H[ l) {k t R)^^H z (p)ds 

+ l -k t P J H[ 1] (k t R)Y {n + [H z {p)] + +n-[H 2 (p)]-} ds 

= [H 8 z (r)] + -±{[H z (r)] + - [ft(r)]_} 

= ^{[H z (r)] + + [HAr)]-) -m(r), reL', (43) 

where L' denotes L less its corners (and endpoints). The same result 
is obtained by letting r tend to ro from the Lj side, as is evident from 
the symmetry in (43). We have made use of the continuity of //Mr) . 
Now n+ = — n_ on L, but we note that the integral equations (41) 
and (43) do not determine H z {r) for r e V, and {[H z (r)]+ — 
[H z (r)]-} for reL', because of the unknown quantity [H z (r)]+ + 
[H z (r)]~ on the right-hand side of (43). Moreover, if L consists of 
segments of a straight line, then R-n = for reL ', p e L' and p ¥>r, 
and the second integral in (43) vanishes. If, in addition, T is empty, 
then (43) reduces to 

[H z (r)] + + [H z (r)]- = 2m(r), reL'. (44) 

This reduction was pointed out by Noble 3 in the case of an infinitesi- 
mally thin strip and by Millar 13 in the case of coplanar strips. In the 
next section, we derive an integral equation which does not degenerate, 
in the case of //waves, for reL'. We remark that the integral equation 
(35) does not degenerate for r e L , and this is presumably because it 
was obtained by setting E\ = —E\ onC. This suggests that we should 
derive an expression for dH 8 z /dn on C, and set it equal to —BH z /dn. 

INTEGRAL EQUATIONS FOR SCATTERING PROBLEMS 417 



IV. INTEGRAL EQUATIONS DERIVED FROM THE GRADIENT OF THE 
SCATTERED FIELD 

We now return to the integral representation (32) for the scattered 
field, and calculate the transverse component of its gradient. If we 
substitute the explicit form (19) of the Green's function into (32), and 
use (20) and (36), we obtain 



•in 



// (i) = 41 I JWiF»<M>5jp+W - HVHkR)^ (p) 



Hence, since 9 

it follows that 

V^"(r) 



d_ 
dR 



[RHFiktR)] = ktRHttHktR), 



ds. (45) 



(46) 



-jAjfflmiO^r^+fr)* 



R ! 



{k t \mHk tR) [l- 



2(Rn)R 

R 3 



^(p)ds 



+ j*< *P(M>55<W* 



Now, since dp/ds = t, and t X i z = — n, we have 



d 



R x U 



-^ + J,(Rxi,)(R.t). 



Also, 



(Rxi 2 )(R-t) = (R-t)[(R-n)t- (R-t)n] 
= (Rn)R-R 2 n. 
Hence, 

n 2(R • n)R d /R X U 



R' 



R 4 



ds\ R- 



(47) 



(48) 



(49) 



(50) 



If we substitute (50) into the second integral in (47), and integrate by 
parts, and combine terms with the help of (49), we obtain 



VW 8 (r) = l -k 2 t J m ) (k t R)n^(p)di 
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4- 



H 



H[ l) (k t R) 
R 



(Rxi 2 )^(p)+R^(p) 
ds on 



ds. (51) 



This expression for the gradient of the scattered field is the 2-dimen- 
sional analog of that derived by Mitzner 4 in 3 dimensions. We give an 
alternate derivation of (51) in the appendix. 

We are interested in calculating the normal derivative of the scat- 
tered field in the limit as r tends to a point r e I" U L', i.e. V,^ 8 (r) • no, 
where no is a unit normal to V U V at ro. It will be seen that the 
limiting value of this quantity may be calculated with the help of (51), 
whereas the second integral in (47) has a singular behavior. However, 
it is not necessary to integrate this second integral by parts completely 
around C, as was done to obtain (51). If we let £ be a segment (or 
segments) of T U L which has r as an interior point, then it suffices to 
integrate by parts over £ . Since the second integral in (47) vanishes 
in the case of E waves, because E z = on the boundary, we now 
consider the case of H waves. 



We define 



= n(r ), roeT', 
110 n+ (r ) = -n_ (r ) , r e U, 



and choose 

n = n+ =* dp/ds = t = n+ X i 2 , on L. 

We also define the tangential component of current density 



J(p) = 



H z (p), peT, 



(52) 



(53) 



(54) 



[H z (p)] + -[H 2 (p)]-, peL. 

Then, from (11) and (47), after an integration by parts, and use of (46), 
(49), and (50), we obtain 



V,/Y 8 2 (r).n = - 



H 

Jr 



HPikR) 



(R.n)(R.no) 



R- 



J(p)ds 



W'iktR) 



n 

R 



2(Rn)R 



R 



n Q J(p)ds 



-k 2 t J m i) (k t R)n-n J(p)ds 
4 k 



4 4 



(Rxi.) dJ 
__ no- (p)ds 
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-L kt \H?Hk t R)?!^.noJ(p)\ 



lo 



(55) 



The contributions from all the endpoints of £ must be included in the 
last term in (55). If £o = T U L, then this last term is zero, and the first 
two integrals in (55) are absent. 

We now consider r — > ro, in a direction which is not tangential to 
I" U U. Since r is an interior point of £ , the first two integrals in 
(55) are well-behaved as r — * ro, as are the contributions from the 
endpoints of £o, represented by the last term in (55). Also, it follows 
from (23) that, for r = ro, there is an integrable singularity in the third 
integral in (55). It remains to consider the fourth integral in (55). As 
depicted in Fig. 3, we take 



r = r + e(cos xto + sin x"o)t 



(56) 



where sin x ** 0- Also, for convenience, we take s = at r = r , so 
that 14 



p = r + sto--K s 2 no + 0(s 3 ) 
for small I s I , where ko is the curvature at ro. Hence, 



(57) 



R = r- p = (ecosx - s)to + (csinx + -Kos 2 )no + 0(s 3 ), (58) 

and 

R 2 = | R 1 2 = ( s - e C os x ) 2 + e 2 sin 2 X + (e s 2 ) + ( s 4 ) . (59) 
From (52), (53), and (58), it follows that 

(Rxi 2 ).no = -R-to = (s- ecosx) +0(s 3 ). (60) 

But 9 




Fig. 3 — Coordinates of a point in the neighborhood of a cross-sectional boundary 
curve. 
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and 

•s 



I 



HPikcR) =-^ + 0[k l R\og(k t R)-\, for &fl« 1, (61) 



(s — e cos x)^ s 



.[(«-€ cob x) 2 + c sin 2 x] 



= 7: log 



(8- ecosx) 2 + c 2 sin 2 x 

.2 . ,2 „• 2 



2 ° (8 + ccosx) +€ sin^x. 



. (62) 



It follows from (59) to (62) that 



lim 



fel 



H^iktR) 



(R x U 



no 



^(p )rf s}=0. 



(63) 



Hence the principal value of the fourth integral in (55) must be taken 
in the limit r — ► To. 

Having shown that the right-hand side of (55) is meaningful in the 
limit r — > ro, we now note that the left-hand side tends to dHl/dno = 
—dHi/dno, since dH z /dn = on the boundary, and hence its value is 
known. Hence the limit of (55) asr— »ro€r'UL' leads to the desired 
integral equation for </(p), as defined in (54). We remark that dJ/ds, 
as well as J, occurs in the integrand. We also remark that, when this 
integral equation has been solved for J(r) for r e T' U L', then (43) 
may be used to calculate [H z (r) ]+ + [H z (r) ]- for r e L' , and hence the 
separate values of [H z (r) ]+ and [H z (r) ]- . We comment that we could 
presumably use the integral equation derived from (55) for ro e L' only, 
and combine it with (41) for r e T', to solve for J(r) for r e V U L'. 

We now consider the case of E waves. Then, from (11) and (47), or 
equivalently (51), 



V,E!(r).no 



H 



R • no dE 2 



\Kt*lf 


R dn ' 


)*)U.O 




R.no f 
R 1 




+ 





ds. (64) 



But, from (58), 

R.no = € sin x + 0(s 2 ), (65) 

and, for 8 > 0, 
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lim 

e*0+ J 



e sin x ds 



~°+ J_, [(* - € cos x) + € sin 2 x] 



= lim 

e-0+ 



tan" 



_i / s — € cos x 
esin x 



= Trsgn(sinx)- (66) 



We first consider r — > r € T'. Then, from (52) and Fig. 3, sin x > 0. 
Hence, from (64), with the help of (59), (61), (65), and (66), we obtain 



4 J Ko an 



H 



H?(kJU) 



Ro* no 
Ro 



BE, 
dn 



w]+[S«]}* 



di?§ . . 1 d£ 2 . . 1 dE z dE z . 

= — - ro)--—(r )=-— -(r )--— (ro), rocT', (67) 

dn 2 dn 2 dn dn 

where Ro = r — p. 

We now consider r — > ro e L'. If the approach is from the plus side 
then, from (52) and Fig. 3, no — n+ (r ) and sin x > 0. Hence, from (64), 
with the help of (59), (61), (65), and (66), it follows that 

Lk t (H[-(k t Ro)^ d -^ip)ds 
4 J ito an 

.r**.,! -i/r*fa.»i +[*w]} 



+ 4 



2[[ dn 



2 l| an 



ds 



(r„)| -|f>> 



}-[§«!■ 



r eL'. (68) 



On the other hand, if the approach to ro is from the minus side, then 
no = — n-(ro) and sin x < 0, and it follows that the left-hand side of 
(68) is equal to 

\ffi\. "I lfTaE,, "1 \dE t . "I 1 

-[-en- iro \ + n^ iro \ + _^r (ro) J_j 
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dn 



dE z 
dn 



(ro) 



]) 



♦[*« 



(69) 



Since [dE z /dn (r )]- = - [dE z /dn (r )]+, for r € U, we again obtain 
(68). 

We remark that when the integral equation (35) has been used to 
determine dE z /dn (r) for r e T' and [dE z /dn (r)] + + [dE z /dn (r)]_ for 
r € L', then (68) may be used to calculate [dE z /dn (r)] + - [dE z / 
dn (r) ]- for r € L', and hence the separate values of [dE z /dn (r) ]+ and 
[dE z /dn (r)]_. This is analogous to the earlier remark concerning the 
use of (43). If L consists of segments of a straight line, then Ro -no = 
for r e L', p e L' and p ¥* r , and the second integral in (68) vanishes. 
If, in addition, T is empty, then (68) reduces to 



\ dE 'i vl 
— (ro) 
dn 


+ 


\ dE *i vl 
— (r ) 
dn 


= 2 


[ dE *i \] 
— (ro) 

dn 



r £ I/. (70) 



This reduction was pointed out by Noble 3 in the case of an inflnitesi- 
mally thin strip. 

It is of interest to note that (35), and the integral equation obtained 
from (55) by letting r -* r e T' U L', are Fredholm equations of the 
first kind. On the other hand, (41) and (67) which hold for r e T', and 
Fredholm equations of the second kind, which is usually preferable 
from the viewpoint of numerical solution. It is somewhat unfortunate, 
in this respect, that the corresponding equations which hold for r e V, 
namely (43) and (68), do not determine [H 2 (r)] + - [H z (r)]- and 
[dE z /dn (r)] + + [dE z /dn (r)]_. 

NOTE ADDED IN PROOF 

In the 3-dimensional case of an open thin shell, the use of a 
degenerate integral equation, analogous to (68), to determine the 
current densities on both sides of the shell, once the total current 
density is known, was pointed out by Stakgold. 16 

APPENDIX 

We here give an alternate derivation of the integral equation (51). 
Let e be a constant vector. Then from (16), it follows that 

(V? + A?)(e.V^')-0, (V? + Jti)(e-W)-0. (71) 

If we apply Green's theorem, as in Section III, but this time to e- V,»p s 
and G, then, with the help of (21) and (27), we obtain 



e.V^'(r) = 



[e-Vty 8 (p)] ^ " G-0- [e.VW(p)] 
dn on 



| Ids, 



(72) 
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where V' ( denotes the transverse component of the gradient with 
respect to the coordinates of p. In a manner analogous to that used in 
Section III, we also obtain 

J ([e.V^'(p)]^-G^[e.V^'(p)]jrfs = 0. (73) 

Hence, with the help of (17), we may rewrite (72) in the form 

e.V„Hr) = J (^Mp)]^ " G ^ [e-VWp)]] da. (74) 
Now 

^■(e.VM-n.VHe.Vty), (75) 

dn 

and 7 

VJ(e.V^) = (e-VO (Vty) = eVfy -Vix(eX Vty) 

= - e*fy -V!x(ex Vty), (76) 

from (16). But 7 

VJ x [G(e x Vty)] = VJG x (e x Vty) + GV;x(ex Vty), (77) 

and, from Stokes' theorem, 15 



I 



n • {V{ X [ G (e X Vty)]} rfs = 0. (78) 

c 

Hence, 

( Gn.[VJ x (e x Vty)]ds = - I n-[V{G x (e x V^)]ds 
Jc Jc 

= - I (n x VfG). (e x Vfflds = I e-[(n x V t G) x V{i//]rfs. (79) 
Jc Jc 



'C JC 

It foUows, from (74) to (76) and (79), that 



-ns" 



e- \V t V(r) - — Vty + tfGnflp) 



Wrjcbj-a 



+ (nxV{G)x Vty \ds \ = 0. (80) 

Since e is an arbitrary (constant) vector, the expression in the curly 
brackets in (80) must vanish. But 

(VlG-nWty + (n X V{G) X Vty 

= (Vty.»)V{G + (n x Vty) x V^G, (81) 
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and n X V t \p = — i 2 bSp/ds on C. Hence, we obtain 



V^ s (r) 



■1 



ki Gn^( P ) + (v;o ^ + (v;g x y ^ 

on ds 



rfs. (82) 



If we substitute for G from (19) and make use of (20) and (36), we 
obtain (51). The representation (82) for the gradient of the scattered 
field is the 2-dimensional analog of that derived by Mitzner 4 in 3 
dimensions. We remark that the above derivation differs from his in 
that we used Green's second identity, whereas he used Green's first 
identity, and consequently some different transformations to reduce 
the result to the form corresponding to (82). Mitzner derived his result 
for the field inside a bounded volume, whereas we are considering 
scattered fields outside bounded cross-sectional curves, and have used 
Green's second identity so that we could make use of the radiation 
condition. 
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